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ABSTRACT 


Propagation  models  used  in  wireless  eommunieation  system  design  play  an  important 
role  in  overall  link  performanee.  Propagation  models  in  a  forested  environment,  in 
particular,  are  especially  valuable  and  complex  due  to  the  randomly  distributed  leaves, 
twigs,  trunks,  and  trees.  This  has  been  an  area  of  interest  due  to  the  operational  needs  of 
military  and  non-military  domains.  Applications  in  both  domains  require  communication 
devices  and  sensors  to  be  operated  in  forested  environments. 

Various  methods  have  been  employed  to  model  propagation  loss.  There  are 
experimental  measurements,  empirical  models,  analytical,  and  computational 
electromagnetic  methods.  Each  method  has  its  applicability  and  limitations.  In  this  thesis, 
investigation  of  a  three-layer  homogenous  medium  model  (air,  forest,  and  ground)  by  a 
ray  tracing  method  was  carried  out.  Both  transmitting  antenna  and  observation  point  are 
within  the  forest  layer. 

The  results  from  the  ray  tracing  model  showed  good  agreement  with  the  available 
measurement  data  up  to  100  MHz.  Further,  better  approximation  of  the  transmission  loss 
was  observed  at  separation  distances  greater  than  1  km.  Values  of  the  effective  electrical 
properties  of  the  forest  played  an  important  role  in  transmission  loss  estimation  within  the 
forest.  After  adjustment,  the  deviation  in  propagation  loss  using  the  ray  tracing  model 
achieved  an  error  of  ±  1  dB. 
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EXECUTIVE  SUMMARY 


Many  applications,  both  in  the  military  and  non-military  domains,  require  the  use  of 
eommunieation  links  within  or  through  foliage.  Examples  of  military  applieations  are  the 
eommunieation  link  between  dismounted  soldiers  and  the  central  node  within  the  forested 
environment,  unmaimed  ground  sensors  (UGS)  deployed  in  a  forested  environment,  and 
downlink  oommunieation  from  an  unmanned  aerial  vehiele  (UAV)  to  troops  loeated 
within  a  forest.  Other  non-military  applieations  are  remote  sensor  networks  for 
agrieulture  monitoring  of  farming  operations  and  seareh  and  reseue  missions  in  a  forested 
environment.  These  real-world  applieations  stimulate  the  need  for  radiowave  propagation 
studies  in  foliage. 

In  an  aetual  forest  environment,  tree  trunks  and  leaves  are  randomly  distributed. 
This  randomness  eauses  attenuation,  scattering,  diffraetion,  and  absorption  of  the  signal 
energy.  Sueh  eharaeteristies  make  the  modeling  of  the  foliage  environment  a  eomplex 
problem.  Moreover,  there  are  various  types  of  woodlands  and  eaeh  has  eharaeteristies 
that  are  different  depending  on  the  elimate.  The  elassifieation  of  different  layers  within 
the  forest  is  shown  in  Figure  1.  This  elassifieation  of  layers  provides  a  simplified  model 
for  analyzing  the  problem  [1].  If  the  variation  within  eaeh  level  as  a  funetion  of  distanee 
is  small  eompared  to  the  wavelength  of  the  propagating  wave,  an  effeetive  homogeneous 
medium  ean  be  defined  for  each  layer.  The  eleetrieal  parameters  that  define  the  media  are 
their  permittivity  and  eonductivity. 


Air  Layer 

omcrQ'cnts  .  •  , 

Canopy  Layer 

. 

Trunk  Layer 

Ground  Layer 

Figure  1.  Fayered  model  of  forest  (after  [1],  [2]). 
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Since  the  1960s,  there  has  been  extensive  research  conducted  to  generate  models 
for  estimating  propagation  loss  [3].  Several  approaehes  have  been  used  to  investigate  and 
model  the  propagation.  The  approaches  inelude  experimental,  empirical,  analytical,  and 
computational  electromagnetic  methods  (CEM).  With  the  advancement  in  teehnology, 
eomputer  proeessors  were  used  to  model  the  propagation  loss.  There  are  several  possible 
eonfigurations  that  occur  when  modeling  the  propagation  loss,  as  depleted  in  Figure  2. 


Air  Air 


Tx  ^ 

Antenna 

j  Forest 

Tx  0 

Antenna _ _ — •• 

X  Forest 

Ground  | 

Ground 

Configuration  1  Configuration  2 
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Antenna  _Air  _ 


Forest 


O 


Ground 


Tx 

Air 

0 

r\lllClllla  1  — ■ 

1  Forest 

Forest  ^ 

Ground 

Configuration  3 


Configuration  4 


Figure  2.  Possible  configurations  to  model  propagation  loss  (after  [4]). 


Eaeh  of  the  approaehes  exhibits  individual  merits  that  contribute  to  a  solution 
addressing  the  problem.  In  the  previous  work  by  Chan  [5],  Computer  Simulation 
Teehnology  (CST)  Microwave  Studio  (MWS)  was  used  for  path  loss  prediction  in 
foliage.  MWS  is  a  three-dimensional  modeling  tool  that  solves  the  integral  form  of 
Maxwell’s  equations  numerically.  It  was  used  to  model  the  eomplex  forest  environment 
as  a  single  layer  dissipative  slab  as  proposed  by  Tamir’s  work  [6].  Simulated  results 
suggest  that  CST  software  ean  be  used  to  provide  a  rough  approximation  of  the 
attenuation;  however,  the  limitation  of  this  technique  is  the  high  eomputational  resources 
needed  and  time  required  for  simulation  of  the  large  slab  (i.e.,  one  beyond  several 
wavelengths). 

In  this  thesis,  an  analytical  method  using  the  ray  traeing  method  is  used  to  model 
the  propagation  with  the  objeetive  of  overeoming  the  shorteomings  of  the  previous  work 
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[5].  The  second  objective  is  to  determine  the  physical  and  electrical  characteristics  of  the 
model  that  correspond  to  the  scenario  of  interest. 

With  the  assumption  that  the  potential  ionoshperic  waves  are  negligible  for  a 
separation  distance  of  less  than  100  km,  two  electromagnetic  (EM)  wave  propagation 
mechanisms  were  identified  within  the  forest  layer  as  depicted  in  Figure  3.  The  two  wave 
components  are  the  lateral  wave  component  (path  TABO)  and  the  Geometrical  Optics 
component  (GO,  direct  and  reflected  waves).  The  formulation  provided  in  [6]  was  used 
for  lateral  component  computation.  Additionally,  modifications  of  the  formulation  used 
to  compute  the  GO  wave  in  [7]  were  made  to  minimize  the  overestimation  of  fields  at 
near  distances.  The  ray  tracing  model,  a  three-layered  model  as  shown  in  Figure  3,  was 
implemented  in  MATFAB.  The  three  layers  defined  in  the  model  are  all  homogeneous. 
Inputs  of  the  geometry  for  the  forest  layer  are  performed  by  the  user.  The  problem 
geometry  setup  is  depicted  in  Figure  4.  Besides  the  geometry  inputs,  the  electrical 
properties  of  each  layer  are  required  from  the  user.  It  was  found  that  only  two  reflections 
from  the  interfaces  need  to  be  considered. 
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Figure  3.  Propagation  mechanism  within  foliage  (after  [6]). 
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Figure  4.  Ray  tracing  model. 


Two  cases  were  analyzed  to  validate  the  ray  tracing  model.  The  objective  of  Case 
I  was  to  compute  the  average  error  between  the  measured  data  extracted  from  [7]  based 
on  the  Sommerfeld  solution  and  the  results  generated  from  the  ray  tracing  model 
(see  Table  1).  The  measured  data  were  taken  by  fixing  the  height  of  the  transmitting 
antenna  to  3.96  m  and  varying  the  height  of  the  receiver  from  7.0  m  to  28.96  m.  The 
separation  distance  D  was  also  kept  constant  at  1 .6  km. 


Table  1 .  Average  error  between  measured  data  from  [7]  and  ray  tracing  model  with 
_ _ adjusted  parameter. _ 


Average  error  compared  between  model  and  measured  data 

H-pol  to  H-pol 

V-pol  to  V-pol 

Frequency 

(MHz) 

Sommerfeld 

solution 

Ray  tracing 
model 

Sommerfeld 

solution 

Ray  tracing 
model 

25 

0.60 

0.82 

0.30 

0.90 

50 

0.80 

0.23 

0.40 

0.48 

100 

0.90 

0.97 

0.40 

0.72 

The  idea  of  using  a  homogenous  medium  to  represent  the  forest  layer  provided 
satisfactory  correlation  to  the  measured  data  (within  1  dB)  up  to  a  frequency  of  100  MHz. 

The  objective  of  Case  II  was  to  compute  the  root  mean  square  (RMS)  error 
between  the  measured  data  from  two  empirical  models  [8]  (the  Jansky  and  Bailey  model, 
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and  Tewari’s  model)  and  the  ray  traeing  model.  A  RMS  error  value  was  introdueed  in  [9] 
to  provide  a  useful  performance  quantity.  In  Case  II,  both  the  heights  of  the  transmitting 
antenna  and  observation  points  were  fixed  at  5.0  m.  The  separation  distance  was  from  1.0 
m  to  4.2  km.  Table  2  is  populated  with  the  resultant  RMS  error  computed  for  each 
deterministic  model  for  distances  greater  than  1 .0  km. 


Table  2. _ RMS  Error  for  measured  data  at  50  MHz  and  D  greater  than  1.0  km. 


Polarization 

Ray  tracing  model, 
^RMS  (dB) 

Tewari’s  empirical 
model, 

^RMS  (dB) 

JB  model, 

^RMS  (dB) 

V-pol 

2.6869 

2.4377 

2.5979 

H-pol 

1.3415 

1.3141 

1.3813 

In  conclusion,  the  ray  tracing  model  is  able  to  provide  good  estimation  of 
propagation  loss  up  to  100  MHz  in  a  forest  environment.  This  model  is  applicable  to 
propagation  loss  estimation  for  cases  where  the  antennas  are  deployed  more  than  1.0  km 
apart.  The  advantage  of  this  model  over  the  previous  model  using  CEM  is  its 
applicability  to  estimate  the  propagation  loss  at  a  distance  beyond  1 .0  km.  Recommended 
future  work  includes  exploring  methods  to  determine  effective  parameters  (permittivity 
and  conductivity)  of  the  forest  layer  and  formulating  models  for  the  other  configurations 
depicted  in  Figure  2.  Lastly,  to  complement  the  empirical  and  analytical  approaches, 
CEM  models  can  be  explored  to  address  propagation  losses  within  the  1 .0  km  range. 


List  of  References 

[1]  L.  W.  Li,  T.  S.  Yeo,  P.  S.  Kooi,  and  M.  S.  Leong,  “Radiowave  propagation  along 
mixed  paths  through  a  four-layered  model  of  rain  forest:  an  analytic  approach,” 
IEEE  Trans,  on  Antennas  Propag.,  vol.  46,  no.  7,  pp.  1096-1  111,  1998. 

[2]  L.  0.  t.  Internet,  “Learn  on  the  Internet,”  [Online].  Accessed  August  2,  2014. 
Available:  http://www.geography.learnonthemternet.co.uk/ 

images/ ecosystems/ crosssection  1  .gif 

[3]  Y.  S.  Meng,  Y.  H.  Lee,  and  B.  C.  Ng,  “Empirical  near  ground  path  loss  modeling 
in  a  forest  at  VHF  and  UHL,”  IEEE  Trans,  on  Antennas  Propag.,  vol.  57,  no.  5, 
pp.  1461-1468,  2009. 


XIX 


[4]  T.  Tamir,  “Radiowave  propagation  along  mixed  paths  in  forest  environments,” 
IEEE  Trans,  on  Antennas  Propag.,  vols.  AP-25,  no.  4,  pp.  471-477,  1977. 

[5]  C.  W.  Chan,  “Investigation  of  propagation  in  foliage  using  simulation 
teehniques,”  M.S.  thesis.  Department  of  Electrieal  and  Computer  Engineering, 
Naval  Postgraduate  Sehool,  Monterey,  CA,  September  2011. 

[6]  T.  Tamir,  “On  radiowave  propagation  in  forest  environments,”  IEEE  Trans,  on 
Antennas  Propag.,  vol.  AP-15,  no.  6,  pp.  806-817,  November  1967. 

[7]  Y.  Ei  and  H.  Eing,  “Numerieal  modeling  and  meehanism  analysis  of  VHE  wave 
propagation  in  forested  environments  using  the  equivalent  slab  model,”  Progress 
in  Electromagn.  Res.,  vol.  91,  pp.  17-34,  2009. 

[8]  R.  K.  Tewari,  S.  Swamp,  and  N.  R.  Manujendra,  “Radiowave  propagation 
through  rain  forests  of  India,”  IEEE  Trans,  on  Antennas  Propag.,  vol.  38,  no.  4, 
pp. 433-449,  1990. 

[9]  M.  O.  Al-Nuaimi  and  R.  Stephens,  “Measurements  and  predietion  model 
optimization  for  signal  attenuation  in  vegetation  media  at  eentimetre  wave 
frequeneies,”  Proc.  in  Microw.,  Antennas  Propag.,  vol.  145,  no.  3,  pp.  201-206, 
1998. 


XX 


ACKNOWLEDGMENTS 


Foremost,  I  would  like  to  thank  Professor  David  C.  Jenn  for  his  guidance  and 
advice.  His  knowledge  and  guidance  was  the  key  to  my  completion  of  this  thesis. 

I  would  like  to  thank  Professor  Phillip  E.  Pace  for  agreeing  to  be  the  second 
reader  of  this  thesis. 

I  would  like  to  thank  my  wife,  Irene  Seah,  for  her  relentless  support  and 
encouragement  during  my  course  of  study.  I  would  also  like  to  thank  my  mother,  Joo 
Khim,  for  her  support  and  encouragement. 

Lastly,  I  would  like  to  thank  Singapore  Technologies  (Electronics),  and  my 
superiors,  Mr.  Chay  Wing  Yew  John  and  Mr.  Ng  Woon  Wei,  for  providing  this 
opportunity  to  pursue  my  master’s  degree  at  the  Naval  Postgraduate  School. 


XXI 


THIS  PAGE  INTENTIONALLY  LEET  BLANK 


I.  INTRODUCTION 


A.  RADIOWAVE  PROPAGATION  IN  A  FORESTED  ENVIRONMENT 

Radiowave  propagation  loss  in  wireless  communication  systems  plays  an 
important  role  in  overall  link  performance.  For  applications  where  communicating 
terminals  are  immersed  in  a  forested  environment,  the  prediction  of  the  propagation  loss 
is  very  challenging.  Tree  trunks,  leaves,  and  undergrowth  in  the  forest  are  randomly 
distributed  and  can  cause  attenuation,  scattering,  diffraction,  and  absorption  of  the  signal 
energy.  This  makes  radiowave  propagation  through  such  environments  a  complex 
problem  to  model.  Modeling  of  this  propagation  loss  has  been  an  area  of  interest  since  the 
1960s  [1]. 

There  are  various  types  of  woodlands  and  each  has  characteristics  that  are 
different  from  one  another  due  to  climate  differences.  A  classification  of  the  different 
layers  within  a  forested  environment  is  depicted  in  Figure  1.  These  layers  provide  a 
simplified  model  for  analyzing  the  problem  [2].  If  the  variation  within  each  level  as  a 
function  of  distance  is  small  compared  to  the  wavelength  of  the  propagating  wave,  then 
an  effective  homogeneous  medium  can  be  defined  for  each  layer. 


Figure  1 .  Layered  model  of  the  rain  forest  (after  [2],  [3]). 
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There  are  many  important  real-world  applieations  that  drive  radiowave 
propagation  studies  in  foliage.  Examples  are  military  operational  scenarios  that  require 
communication  systems  to  be  deployed  in  forested  environments,  including  unmanned 
ground  sensors  (UGS)  and  downlink  communication  from  an  unmanned  aerial  vehicle 
(UAV)  to  troops  located  within  forest.  Lastly,  there  is  deployment  of  vehicular 
communication  nodes  within  forested  environments  to  relay  information  to  dismounted 
soldiers  [4]. 

In  the  context  of  non-military  operations,  there  are  search  and  rescue  missions  in 
forests  [5]  and  also  wireless  sensor  networks  for  agricultural  applications  [6].  These  are 
just  some  examples  of  communication  links  that  require  consideration  of  the  losses  due  to 
the  scattering,  diffraction,  and  attenuation  within  the  forested  environment. 

B,  BACKGROUND 

Several  approaches  have  been  used  to  investigate  and  model  the  propagation  loss 
characteristics  of  a  complex  forested  environment.  One  of  these  approaches  is 
experimental  measurement.  Numerous  radiowave  measurement  studies  have  been 
conducted  in  the  tropical  rainforest  located  in  Northern  India,  where  the  measured  results 
were  taken  in  the  frequency  range  from  50  MHz  to  800  MHz  [7].  Other  measurements 
were  taken  from  central  and  southern  parts  of  Thailand  [8]  at  frequency  ranges  from  2 
MHz  to  400  MHz.  The  limitation  to  these  measurements  is  that  they  cannot  easily  be 
extended  to  other  types  of  forest. 

Empirical  models  were  also  introduced  for  predicting  the  propagation  loss  from 
various  sources.  They  include  the  International  Telecommunication  Union 
Recommendations  (ITU-R)  [9],  Modified  Exponential  Decay  (MED)  model  by 
Weissberger  [10],  and  COST  235  model  [11].  There  are  many  other  models  in  addition  to 
these.  The  models  are  useful  in  providing  a  rough  estimate  of  radiowave  attenuation 
through  vegetation. 

Analytical  methods  have  been  used  to  determine  the  characteristics  of  radiowave 
propagation  in  a  forested  environment.  Tamir  [12]  proposed  that  the  complex  forested 
environment  can  be  modeled  as  a  dissipative  slab  provided  the  wavelength  at  the 
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operating  frequency  is  much  smaller  than  the  irregularities  of  the  forested  environment. 
This  model,  together  with  the  theory  on  wave  propagation  in  a  stratified  medium,  led  to 
the  finding  that  the  “tree  top”  wave  propagation  is  the  dominant  mode  compared  to  the 
direct  and  multipath  waves  within  the  forest.  Leveraging  this  concept,  there  were 
subsequent  works  done  using  the  same  idea  of  treating  the  forested  layers  as 
homogenous,  isotropic  media  [13],  [14]. 

With  the  advancement  in  computer  processors,  the  use  of  computational 
electromagnetic  methods  (CEM)  has  become  an  attractive  tool  to  solve  more  complex 
real-world  electromagnetic  (EM)  problems.  Eeveraging  on  this,  the  previous  work  by 
Chan  [13]  investigated  simulation  techniques  for  path  loss  prediction  in  foliage  using 
Computer  Simulation  Technology  (CST)  Microwave  Studio  (MWS).  MWS  is  a  three- 
dimensional  modeling  tool  that  solves  the  integral  form  of  Maxwell’s  equations 
numerically.  It  was  used  to  model  the  complex  forest  environment  as  a  single  layer 
dissipative  slab  as  proposed  by  Tamir’s  work.  Simulated  results  suggest  that  the  model  of 
the  slab  using  the  CST  software  can  be  used  to  provide  a  rough  approximation  of  the 
attenuation.  The  limitation  of  this  technique  is  the  high  computational  resources  needed 
and  time  required  for  simulation  of  large  slab  (i.e.,  one  beyond  several  wavelengths). 

Other  approaches  have  been  used  to  model  the  radiowave  characteristics  in  a 
forested  environment.  The  frequency  range  of  the  early  works  is  the  very  high  frequency 
(VHE)  band  (30  to  300  MHz)  where  the  wavelength  is  relatively  long  (1.0  m  to  10.0  m). 
With  the  congestion  of  spectrum  usage  in  the  VHE  band  and  the  introduction  of  digital 
communication  systems  into  the  communication  community,  there  has  been  a  shift  of 
frequency  usage  to  the  ultra-high  frequency  (UHE)  band  [15].  Operating  at  higher 
frequencies  in  the  forest  environment  poses  different  challenges,  as  the  wavelengths  are 
shorter  for  higher  frequencies.  The  work  done  by  Ei  [14]  is  of  particular  interest.  The  use 
of  ray  tracing  methods  suggests  that  the  numerical  results  for  the  transmission  loss  and 
the  time  delay  versus  range  exhibit  good  agreement  with  the  numerical  Sommerfeld 
integral  solutions.  Eurther  modeling  of  the  propagation  using  the  ray  tracing  method  is 
presented  in  this  thesis. 
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c. 


THESIS  OBJECTIVE 


In  the  previous  study  by  Chan  [13],  MWS  was  used  to  model  a  homogeneous 
lossy  dielectric  slab  as  a  representation  of  the  forest  medium.  Two  important  parameters 
(conductivity  and  permittivity)  used  in  modeling  the  slab  that  determine  the  forest 
electrical  characteristics  are  discussed  in  this  thesis.  A  rough  approximation  for 
radiowave  propagation  at  short  distances  within  the  forest  medium  was  achieved.  One  of 
the  limitations  in  using  CEM  to  model  the  radiowave  propagation  loss  is  the  long 
processing  time  and  large  amount  of  computer  memory  required.  One  objective  of  this 
thesis  is  to  explore  an  alternative  approach  that  removes  this  limitation. 

A  deterministic  model  based  on  ray  tracing  methods  is  investigated  in  this 
research.  This  model  is  used  to  generate  the  lateral  waves  as  well  as  the  geometric  optic 
(GO)  waves  which  are  summed  at  the  observation  point,  as  illustrated  in  Figure  2. 
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Figure  2.  Illustration  of  propagation  mechanisms  (after  [12]). 


Simulation  results  from  the  code  generated  in  MATFAB  software  are  used  to 
compare  with  measured  data  extracted  from  [7]  and  [14].  This  is  to  re-establish  the 
finding  in  [12]  that  a  homogeneous  slab  is  applicable  to  compute  propagation  loss  up  to 
100  MHz  within  the  VHF  band. 
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D,  THESIS  OUTLINE 

There  are  five  chapters  in  this  thesis  and  they  are  organized  in  the  following 
sequence. 

In  Chapter  I,  the  background  and  objectives  of  the  thesis  was  presented. 

In  Chapter  II,  the  EM  wave  propagation  mechanisms  within  the  forest  layer,  as 
well  as  the  formulations  for  the  different  mechanisms  used  to  describe  the  wave 
propagation,  are  discussed. 

In  Chapter  III,  we  document  the  list  of  assumptions  and  considerations  for  the  ray 
tracing  model.  Formulation  of  the  model  are  described  in  detail. 

In  Chapter  IV,  comparisons  of  the  simulated  results  are  made  with  measured  data 
and  some  observations  are  offered  based  on  the  results. 

Finally,  in  Chapter  V,  a  summary  and  conclusions  with  recommendations  for 
future  work  are  presented. 
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II.  METHODOLOGY  USED  FOR  MODELING 


In  this  chapter,  the  various  analytical  methods  employed  to  study  the  propagation 
loss  for  communication  links  within  the  forest  in  the  VHP  and  UHF  bands  are 
summarized.  A  set  of  possible  eonfigurations  for  communieation  systems  operating 
within  or  outside  a  forest  environment  is  defined.  A  summary  of  EM  plane  wave  theory  is 
provided  to  illustrate  the  effeet  of  various  parameters. 

Besides  reviewing  the  analytical  methods,  the  propagation  mechanisms  that  are 
present  within  the  forest  environment  are  described.  A  list  of  empirieal  models  used  in 
the  study  of  propagation  loss  from  foliage  are  presented,  to  be  used  as  part  of  the  tools  to 
validate  our  model 

A,  ANALYTICAL  METHODS 

Analytical  approaches  have  been  used  in  modeling  propagation  loss. 
Deterministie  models  use  different  solution  techniques  to  eompute  the  propagation  loss 
for  the  various  eonfigurations  depieted  in  Figure  3  [16]-[19]. 
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Figure  3.  Various  configurations  for  deterministic  modeling  (after  [18]). 
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Most  of  the  analytical  modeling  uses  an  infinitesimal  dipole  as  the  transmitting 
antenna  or  for  both  transmitting  and  receiving.  The  point  O  labeled  in  Figure  3  is  either 
the  location  of  a  receiving  antenna  or  simply  a  point  in  space  for  computation  of  electric 
field  strength.  The  calculation  of  received  power  requires  that  an  antenna  be  specified; 
however,  the  calculation  of  field  strength  does  not. 

The  initial  work  by  Tamir  suggested  that  the  upper  limit  of  the  operating 
frequency  for  the  foliage  block  model  is  100  MHz  [12].  In  his  later  work  with  Dence 
[20],  this  upper  limit  was  extended  to  200  MHz.  This  extension  of  upper  frequency  limit 
was  obtained  by  including  the  ground’s  effect  into  the  initial  closed  form  lateral  wave 
equation  derived  in  [12].  A  different  method  was  reported  by  Palud  [21]  in  2004,  and  he 
compared  his  method  used  with  the  results  reported  in  [18].  Palud  used  the  full-wave 
Parabolic  Equation  (PE)  algorithm  in  treating  the  forest  as  a  lossy  medium  characterized 
by  its  permittivity  and  conductivity.  The  model  was  able  to  predict  the  propagation  loss 
in  either  Configuration  1  or  Configuration  2  depicted  in  Eigure  3.  The  results  for 
Configuration  1  showed  noticeable  asymmetry  in  the  uplink  and  downlink  and  did  not 
match  to  the  results  that  were  reported  in  [18];  however,  results  based  on  Configuration  2 
showed  good  agreement.  Eiao  and  Saraband!  [16]  investigated  the  problem  by  using  a 
surface  field  integration  technique.  Configuration  1  was  studied,  and  it  was  found  that  at 
low  observation  heights,  the  solution  underestimates  the  received  field  at  the  observation 
point  compared  to  [18]  due  to  the  assumption  made  by  Tamir  that  all  energy  arriving  at 
the  canopy  edge  is  assumed  to  be  coupled  entirely  towards  the  observation  point  for  both 
polarizations. 

Given  the  same  idea  of  a  homogeneous  medium  to  represent  the  forest  layer,  there 
are  various  deterministic  models  proposed  for  planar  stratified  media.  A  three-layered 
model  using  Green’s  dyadic  function  to  predict  propagation  loss  was  reported  in  [22]. 
The  theoretical  results  were  in  good  agreement  with  the  available  experimental  data  at 
6  MHz.  Within  the  same  time  frame,  Cavalante  used  the  same  method  and  modeled  the 
propagation  loss  using  four  homogenous  layers  [23].  This  four-layer  model  takes  into 
account  the  vertical  non-homogeneities  of  the  forests  in  order  to  extend  the  prediction 
accuracy  as  high  as  250  MHz.  Subsequently,  Seker  and  Scheider  [15]  used  the  Hertz 
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potential  method  applied  to  anisotropie  stratified  media  to  prediet  the  propagation  loss  in 
a  forest  environment.  The  range  of  frequeneies  was  from  20  MHz  to  2,000  MHz  due  to 
the  high  level  of  interest  in  several  wideband  digital  radio  systems  [24],  However,  results 
from  both  [15]  and  [24]  were  not  eompared  to  any  measured  data  at  the  time. 

Extensive  analytieal  study  on  the  propagation  loss  at  both  VHP  and  UHF  were 
eondueted  by  Li  [1].  This  extensive  study  of  radiowave  propagation  in  a  four-layered 
model  uses  Green’s  dyadie  funetion  to  solve  the  problem.  Li  eoneluded  that  the  lateral 
waves  at  the  forest-air  boundary  are  the  dominant  waves  eompared  to  the  other  lateral 
waves  propagating  along  the  air-eanopy,  trunk-eanopy  and  ground-trunk  boundaries.  The 
high  attenuation  of  lateral  waves  along  the  ground-trunk  boundary  is  due  to  the  high 
eonduetivity  of  the  ground  layer  [19]. 

A  summary  of  the  various  solution  approaches  and  the  applicable  configurations 
defined  in  Figure  3  are  tabulated  in  Table  1. 


Table  1 .  Summary  table  for  the  various  analytic  methods  used  to  model 

propagation  loss. 


Method 

Configuration  with 
reference  to  Figure  3 

Number  of 
Layers 

Planar  Stratified 
Mediums 

Parabolic  Equation 
Algorithm  [21] 

1,2 

2 

Homogeneous 

Ray  Tracing  [18] 

2 

2 

Homogeneous 

Sommerfeld  Integral 
Solution  [14] 

2 

3 

Anisotropic 

Dyadic  Green’s 
Function  [19]  ,  [22] 
[23]  and  [25] 

2 

3,4 

Homogeneous 

Surface  Field  Integral 
Technique  [16] 

1 

2 

Homogeneous 

Hertz  Potential 
Method  [15]  and  [24] 

2 

4 

Anisotropic 
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We  use  a  ray  tracing  method  to  model  a  three-layer  planar  stratified  medium  as 
shown  in  Figure  2,  and  the  details  of  the  model  are  covered  in  Chapter  III.  For  our 
analytical  approach,  we  require  measured  or  estimated  electrical  data  as  inputs  for 
properties  of  the  forest  medium.  A  review  of  the  various  approaches  used  to  obtain  the 
electrical  properties  of  the  forest  medium  was  performed  and  the  values  reported  in  [26]. 
This  is  discussed  in  Section  C  of  this  chapter. 

B.  EM  PLANE  WAVE  PROPAGATION 

EM  waves  consist  of  two  coupled,  time-varying  vector  fields,  E  (electric)  and  H 

(magnetic).  The  time  dependence  is  assumed  in  our  analysis  and  suppressed 

throughout  this  thesis.  These  fields  are  governed  by  physical  laws  described  by 
Maxwell’s  equations.  As  an  illustration  of  the  fields  in  space,  let  the  direction  of 
propagation  of  the  wave  be  in  the  positive  z -direction.  The  E  fields  and  H  fields  are 
pointing  in  the  positive  x-  and  positive  y -directions,  respectively,  as  shown  in  Figure  4. 


Plane  waves  are  assumed  in  our  discussion  of  wave  propagation.  The  wave  fronts 
are  infinite  parallel  planes  of  constant  amplitude  normal  to  the  phase  velocity  vector. 
Plane  wave  approximations  are  adequate  when  the  extent  of  the  boundaries  between  the 


10 


layers  are  large  eompared  to  wavelength.  Then  the  refleetion  and  refraetion  of  a  spherical 
wave  can  be  considered  locally  planar. 

In  lossless  media,  Maxwell’s  equations  in  phasor  form  relate  the  .Eand  H  fields 

as  [27]; 

VxE  =  -jo)juH  (1) 

VxH  =  -jcosE  (2) 

where  co  =  In  f  is  the  angular  frequency  in  rad/s,  //  is  the  permeability  of  medium,  e  is 
the  permittivity  of  the  medium,  and  j  =  Vm  . 

The  wave  equations  are  derived  from  Eq.  (1)  and  Eq.  (2),  resulting  in  the  final 
form  expressed  as  [27]; 

+  o?/ueE  =  0  (3) 

+  ap-jusH  =  0  .  (4) 

Now  the  wavenumber  k  of  the  medium  is  expressed  as; 

k  -  (5) 

where  values  of  and  s  are  real  in  a  lossless  medium. 


Eor  plane  waves  in  a  lossy  medium.  Maxwell’s  equations  can  be  written  as  [27]; 

V  xE  =  -jcoiuH  (6) 

WxH  =  jcosE  +  aE  (7) 

where  a  is  the  conductivity  of  the  medium,  s  -  s^Sq  where  s  is  the  permittivity  and  Sy. 


is  the  relative  permittivity.  The  complex  permittivity  is  defined  as  [12]; 


Sc  -  SQSy  -  SQ 


^r-J- 


COSq 


(8) 


where  Sy  is  the  complex  relative  permittivity  of  the  lossy  medium.  The  resulting  wave 
equation  for  E  is  now  expressed  as; 

w'^E  +  Qp  /us 


1-7  — 
cos 


E  =  Q  . 


(9) 
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The  difference  between  Eq.  (3)  and  Eq.  (9)  is  that  the  complex  permittivity  introduced 

by  the  lossy  medium  makes  the  wavenumber  a  complex  number.  This  complex  number  is 
defined  as  the  propagation  constant  y  expressed  as: 


y  =  a  +  jP  =  j(Uyf/j£  I  -j  — 

V  coe 

where  a  is  the  attenuation  constant  and  is  the  phase  constant. 


(10) 


With  the  assumption  that  the  wave  is  propagating  in  the  positive  z -direction, 
Eq.  (9)  can  be  reduced  to: 


dz^ 


(11) 


and  the  solution  for  E'^(z)  is: 


E^{z)^E^e~^^  +E~e^^.  (12) 

The  exponential  factor  for  travelling  waves  in  the  positive  z -direction  is  expressed  as: 

e-r  (13) 

The  ratio  between  the  E  and  H  fields  determines  the  intrinsic  impedance. 
Computed  value  for  impedance  in  free  space  is  rjQ  =  ^/uq  /  Sq  =  (120;t)  Q.  .  The 
impedance  of  the  forest  is  expressed  as  [27]: 

120;r 


/  / -  izu;t 

^/=Vv^/'-c  "1— 

rf'^c 


(14) 


where  is  relative  complex  permittivity  of  the  forest  layer  expressed  as  [12] 


where  Sj-  is  the  relative  permittivity  of  the  forest  layer,  ct/  is  the  conductivity  of  the 
forest  layer  and  A  is  the  wavelength  in  meters. 


Until  this  point,  we  have  covered  the  theory  required  to  compute  E  within  a 
medium  with  a  given  dielectric  constant  and  conductivity.  To  address  the  waves  reflected 
off  the  boundaries  of  the  medium,  there  is  another  parameter  to  be  considered.  The 

computation  of  the  reflection  coefficient  for  waves  reflected  off  of  a  plane  interface 
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between  two  dielectric  regions  is  depicted  in  Figure  5.  Based  on  the  well-known  Snell’s 
laws  of  reflection  and  refraction  [27]; 

Oi-0,  (16) 

sin  6i  -  rj2  sin  9^  (17) 

where  Of  is  the  incident  wave  angle  from  the  normal  to  the  plane  interface,  9j,  is  the 

reflected  wave  angle,  Of  is  the  transmitted  wave  angle  or  refracted  angle.  Lastly,  //j  and 

72  are  the  impedances  for  medium  1  and  medium  2,  respectively. 


A  A 

X  X 


Figure  5.  Polarization  for  electric  field  vector  (after  [27]). 


The  polarization  for  the  EM  wave  is  with  reference  to  the  plane  of  incidence 
{x-z  plane)  as  shown  in  Figure  5.  For  vertical  polarization  (V-pol),  the  vector  Ey  is 

parallel  to  the  x-z  plane.  The  component  of  Eh  perpendicular  to  the  x-z  plane  is 
defined  as  horizontal  polarization  (H-pol).  These  terms  originate  from  low  grazing  angles 
—>90°) .  The  parallel  component  becomes  vertical  with  respect  to  the  interface.  The 
perpendicular  component  is  horizontal. 

The  reflection  coefficient  for  vertical  polarization  is  expressed  as  [27]: 

7l  COS  6^^  +  72  COS  6]' 

and  the  reflection  coefficient  for  horizontal  polarization  is  expressed  as  [27]; 
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(19) 


^  _  72  cos  6'- -7/1  cos  6'^ 

^  1  -  ^  if  - - TT^ - IT 

72  cos  +  7i  cos 

where  %  and  72  are  the  wave  impedanees  for  medium  1  and  medium  2,  respectively, 
and  Of  is  related  to  Of  by  Eq.  (17). 

C.  ELECTRICAL  PROPERTIES  OF  THE  MEDIA 

A  summary  of  the  various  approaehes  used  in  determining  the  eleetrieal 
properties  of  the  forest  is  deseribed  in  [26].  Open-wire  transmission  lines  (OWLs)  were 
used  as  probes  to  make  actual  site  measurements  of  the  eleetrieal  properties  of  a  forest  by 
Parker  and  Hagn  [28]  in  1966.  The  study  using  the  OWL  method  was  restricted  to 
75  MHz.  Both  Parker  and  Hagn  also  eondueted  a  feasibility  study  on  using  the  parallel 
plate  eapaeitor  teehnique  and  eavity  teehnique  to  estimate  the  eleetrieal  properties  of  the 
forest  in  the  same  study  published  in  1966.  The  parallel  plate  eapaeitor  teehnique  was 
found  to  have  many  limitations,  and  the  eavity  teehnique  was  not  suitable  [26].  Another 
approaeh  introduced  to  determine  these  eleetrieal  properties  is  the  inverse  method.  This 
method  utilizes  the  data  eolleeted  from  measurements  sueh  as  antenna  patterns  within 
vegetation,  transmission  loss  over  distanee,  wave  impedance,  or  attenuation  variations  to 
estimate  the  values  of  relative  permittivity  and  eonduetivity  ct/  of  the  forest  slab. 

The  eleetrieal  properties  of  the  forest  medium  obtained  via  the  inverse  method  from  [26] 
are  tabulated  in  Table  2.  Eleetrieal  properties  from  eaeh  site  are  different,  as  expeeted,  as 
well  as  frequeney  dependent. 
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Table  2.  Electrical  properties  of  forest  medium  (after  [26]). 


Measurement 

Site 

Source 

Frequency 

(MHz) 

Relative 

permittivity 

^fr 

Conductivity 

(mS/m) 

Attenuation 
constant  cry 

(Np/m) 

India 

(Northern 

India) 

Tewari 

(1982) 

50 

1.065 

0.135 

0.0246 

200 

1.055 

0.145 

0.0266 

500 

1.040 

0.160 

0.0296 

800 

1.040 

0.160 

0.0296 

Pak  Chong 

(Northern 

Thailand) 

Jansky  and 

Bailey 

(1966-73) 

2 

1.010 

0.045 

0.0083 

400 

1.010 

0.045 

0.0084 

Satun 

(Southern 

Thailand) 

Jansky  and 

Bailey 

(1966-73) 

2 

1.010 

0.035 

0.0065 

400 

1.010 

0.035 

0.0066 

D,  ELECTROMAGNETIC  WAVE  PROPAGATION  MECHANISM  FOR 
FOREST  LAYER 

The  two  main  propagation  mechanisms  used  to  characterize  the  communication 
link  within  the  forest  layer  are  lateral  waves  and  Geometrical  Optics  (GO)  waves  as 
described  in  [14],  Lateral  waves  are  found  propagating  along  the  air-forest  boundary 
(path  AB)  depicted  in  Figure  2.  This  lateral  wave  is  generated  by  the  refraction  of  the  GO 
wave  incident  at  a  critical  angle  and  propagates  along  the  path  AB.  As  it  propagates 

along  the  forest-air  boundary,  it  radiates  back  into  the  forest  layer  at  the  same  critical 
angle. 


This  critical  angle  is  dependent  on  the  relative  permittivity  of  the  forest 


defined  as  [12]; 


=  sin  ^ 


^  1  ^ 


y^frj 


(20) 


The  real  part  of  the  complex  permittivity  of  the  forest  layer  is  used  in  Eq.  (20).  To  obtain 
the  lateral  wave  component  travelling  along  the  forest-air  boundary,  the  GO  wave 
traveling  towards  the  forest-air  boundary  at  this  critical  angle  assumes  the  imaginary  term 
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of  the  complex  permittivity  is  very  much  smaller  compared  to  the  magnitude  of  the 
complex  permittivity  [29]. 

Multiple  reflected  GO  waves  are  generated  within  the  forest  layer  either  via 
reflections  off  the  forest-air  boundary  or  the  ground-forest  boundary.  The  effect  of  the 
reflected  GO  waves  with  more  than  two  reflections  off  the  same  boundaries  has  little 
effect  on  the  overall  electric  field.  There  are  two  main  contributing  factors  to  this  claim. 
The  first  is  the  high  attenuation.  Let  us  take  the  India  forest  example  from  Table  2.  For 
50  MHz  a  direct  line  -  of  -  sight  communication  link  with  a  distance  of  100  meters  apart, 
the  attenuation  constant  is  about  0.0246  Np/m.  At  this  distance  (100  m),  the  total  loss  is 
computed  by  multiplying  the  attenuation  constant  value  with  a  conversion  factor  of 
8.686.  This  results  in  the  total  attenuation  of  21  dB. 

The  second  contributing  factor  is  the  reflection  coefficient.  There  is  multiplication 
of  the  reflection  coefficient  factor  in  addition  to  the  attenuation.  Each  reflection 
coefficient  multiplication  reduces  the  field.  The  field  is  further  attenuated  by  the 
additional  path  traveled  in  the  lossy  foliage  to  the  next  reflection  point;  thus,  we  are  only 
taking  into  account  waves  that  have  reflections  off  the  boundaries  at  most  once. 

E,  EMPIRICAL  MODELS 

In  this  section,  seven  models  are  discussed;  the  Weissberger  model.  International 
Telecommunication  Union  -  Recommended  (ITU-R)  model,  fitted  ITU-R  model,  COST 
235  model,  lateral  ITU-R  model,  Jansky  and  Bailey  model,  and  Tewari’s  empirical 
model.  The  total  propagation  loss  through  the  vegetation  for  some  of  the  models 
(Weissberger,  ITU-R,  Fitted  ITU-R,  COST235,  and  Lateral  ITU-R)  requires  the  addition 
of  the  plane  earth  attenuation  PE  defined  as: 

P£(dB)  =  40  logio  ( L>)  -  20  logl  0  ( ^ )  “  20  logj  o  ( /?2 )  (21) 

where  ky  is  the  height  of  transmitting  antenna  in  meters  and  ^2  is  the  height  of 
observation  point  in  meters.  The  separation  distance  D  is  assumed  to  be  very  large 
compared  to  hi+h2  [30]. 
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(1)  Weissberger  Model 

The  Weissberger  model  was  introduced  in  1982  to  provide  a  better  prediction 
from  the  exponential  decay  (EXD)  model  [10].  In  the  report,  it  was  called  the  modified 
exponential  decay  (MED)  model  and  is  applicable  for  an  operational  scenario  where  both 
transmitting  and  receiving  antennas  are  immersed  in  temperate-latitude  forests  filled  with 
dense,  dry,  in-leaf  trees.  The  model  predicts  the  foliage  loss  excluding  the  plane  earth 
losses.  The  frequency  can  be  from  230  MHz  to  95  GHz.  Therefore,  addition  of  the  plane 
earth  model  losses  is  required. 


The  empirical  model  is  expressed  as: 


14  m  <D  <400m 
Om  <  D  <  14m 


(22) 


where  L^r  is  the  vegetation  loss  in  dB,  is  the  frequency  in  GHz  and  D  is  the 

separation  distance  in  meters  between  the  transmitting  antenna  and  observation  point. 


(2)  International  Telecommunication  Union  -  Recommendation  (ITU-R) 

Model 


The  ITU-R  model  [9]  was  developed  from  measured  data  in  the  UHE  band  in 
1986.  In  measurement  setup,  the  transmitting  antenna,  and  observation  point  are 
separated  by  a  small  grove  of  trees  with  depth  less  than  400  m.  The  frequency  can  be 
from  200  MHz  to  95  GHz.  The  majority  of  the  signal  propagates  through  the  grove  of 
trees,  and  the  empirical  model  is  expressed  as: 

L,TU-R(dB)  =  0.2f'>^^D'‘<’  (23) 

where  Litu-R  vegetation  loss  in  dB  and  the  frequency  /  is  in  MHz. 


(3)  Fitted  ITU-R  Model 

This  model  was  proposed  by  Al-Nuaimi  and  Stephens  [31]  from  optimization  of 
the  ITU-R  model  by  measurements  made  at  11.2  GHz  and  20  GHz  at  two  foliation  states 
(in-leaf  and  out  of  leaf).  The  expression  for  the  empirical  model  is: 
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^FITU-R 


out  of  leaf 
in  leaf 


where  Lpij’jj_ji  is  the  vegetation  loss  in  dB  and  the  frequency  /  is  in  MHz. 


(24) 


(4)  COST  235  Model 


The  COST  235  model  [11]  is  proposed  based  on  measured  data  between  9.6  GHz 
and  57.6  GHz  through  a  small  grove  of  trees  less  than  200  m.  The  empirical  model  is 
expressed  as: 


Lc0STidB)-\ 


26.6/“®-^Z)®'^ 

15.6/“®'®®^/)®'^^ 


out  of  leaf 
in  leaf 


(25) 


where  LpoST  the  vegetation  loss  in  dB  and  the  frequency  /  is  in  MHz. 


(5)  Lateral  ITU-R  Model 

The  lateral  ITU-R  (LITU-R)  model  was  proposed  by  Meng  [32]  to  analyze  the 
effect  of  foliage  and  the  ground’s  effect  on  the  propagation  loss  over  a  tropical  plantation 
at  both  VHP  and  UHF  bands.  The  modified  ITU-R  model  is  express  as: 

LuTU-R(‘IB)s0.4sf''^D'>'^  (26) 

where  Liipjj_p  is  the  vegetation  loss  in  dB  and  the  frequency  /  is  in  MHz. 


(6)  Jansky  and  Bailey  Model 

The  Jansky  and  Bailey  (JB)  model  was  developed  based  on  the  measurement  data 
from  a  tropical  forest  in  Pak  Chong,  Thailand  [10].  The  model  applies  to  Configuration  2 
in  Figure  3  where  both  the  transmitting  antenna  and  observation  point  are  immersed 
within  the  forest  layer.  The  frequency  can  be  from  25  MHz  to  400  MHz.  It  can  be 
adapted  to  either  polarization  with  the  heights  for  both  transmitting  and  receiving 
antennas  between  2.0  and  7.0  m.  The  range  of  the  separation  distance  D  is  defined  from 
8.0  m  to  1.6  km. 
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The  loss  expression  for  the  JB  model  is  [10]: 


Ljb  (dB)  =  36.57  +  20  logjo  /  -  20  logjo 


-1609a  [rD. 


D 


JB^mile  Q 
+ 


JB 


mile 


r)2 

^mile 


(27) 


where  is  the  separation  distanee  between  the  transmitting  and  reeeiving  antennas 

in  miles  and  Ljg  is  the  total  transmission  loss  in  dB  and  the  frequeney  /  is  in  MHz. 
^JB  and  Bjb  are  eonstants  obtained  empirieally  from  the  measured  data.  These 
eonstants  are  listed  in  Table  3. 


Table  3.  Constants  for  the  JB  Model  (after  [10]). 


Frequency 

(MHz) 

Polarization 

^JB 

(dB/m) 

^JB 

Bjb 

25 

V-pol 

0.0 

0.0 

0.002120 

50 

V-pol 

0.0 

0.0 

0.001060 

100 

V-pol 

0.045 

0.615 

0.000529 

250 

V-pol 

0.050 

0.759 

0.000443 

400 

V-pol 

0.055 

1.020 

0.000523 

25 

H-pol 

0.0 

0.0 

0.004240 

50 

H-pol 

0.0 

0.0 

0.004240 

100 

H-pol 

0.020 

0.472 

0.005510 

250 

H-pol 

0.025 

0.774 

0.000588 

400 

H-pol 

0.035 

1.110 

0.000598 

(7)  Tewari’s  Empirical  Model 

Tewari’s  empirieal  model  was  developed  using  the  measured  data  made  in  a 
reserve  forest  near  Dehradun,  Northern  India  [7].  The  model  addresses  the  exponential 
decay  characteristics  for  forest  depths  less  than  400  m.  It  was  suggested  that  beyond  the 
nominal  400  m,  received  signals  from  the  transmitting  antenna  were  partially  via  the 
lateral  mode.  Similar  to  that  of  the  JB  model,  there  is  also  a  table  of  constants,  shown  in 
Table  4,  applied  for  V-pol  and  H-pol.  The  applicable  range  of  this  model  is  from 
40  m  to  4  km.  The  frequency  range  is  from  50  MHz  to  800  MHz. 
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Table  4.  Constants  for  the  Tewari’s  empirieal  model  (after  [7]). 


Frequency 

(MHz) 

Polarization 

Uj 

(dB/m) 

Ay 

By 

50 

V-pol 

0.0 

0.0 

1.9170 

200 

V-pol 

0.0125 

0.4989 

1.8358 

500 

V-pol 

0.0135 

0.3658 

0.9040 

800 

V-pol 

0.0140 

0.2661 

0.5331 

50 

H-pol 

0.0 

0.0 

7.3670 

200 

H-pol 

0.0110 

0.8201 

5.0450 

500 

H-pol 

0.0138 

0.6571 

1.4304 

800 

H-pol 

0.0152 

0.4491 

0.6291 

The  loss  expression  for  Tewari’s  empirieal  model  is  [7]; 


Lj  {dB)  =  -27.57  +  20  logj  q  /  -  20  logig 


D 


By 

D 


where  Lj  is  the  total  transmission  loss  in  dB  and  the  frequeney  /  is  in  MHz. 


(28) 


In  summary,  we  have  deseribed  various  eonfigurations  for  deterministie  modeling 
of  propagation  loss  due  to  the  presenee  of  vegetation.  A  summary  of  the  various 
approaehes  with  referenee  to  [1]  was  presented.  The  various  analytieal  models  presented 
were  foeused  on  the  study  of  lateral  waves  and  the  waves  generated  in  the  VHP  and  UHF 
bands.  Multiple  planar  media  were  utilized  to  eharaeterize  the  different  layers  of  the 
forest  like  the  trunk  and  eanopy  layers.  As  the  frequeney  of  interest  is  inereased  to 
millimeter  waves,  the  lateral  wave  over  large  forests  ean  no  longer  exist.  This  leads  to  the 
use  of  CEM  to  provide  solutions  to  the  intraetable  mathematieal  formulations  but  also 
demands  more  eomputational  resourees. 

The  ray  traeing  method  used  to  prediet  the  propagation  loss  within  a  forest  was 
introdueed.  We  shall  make  use  of  this  knowledge  to  formulate  the  model  in  Chapter  III. 
Different  approaehes  were  used  to  determine  the  eleetrieal  properties  of  the  forest  layer. 
These  eleetrieal  properties  are  site  speeifie,  and  eleetrieal  properties  for  two  sites  were 
given  in  Table  2.  This  table  is  eonsulted  again  in  our  analysis  in  Chapter  IV.  Finally,  we 
have  presented  the  empirieal  models  used  to  prediet  the  propagation  loss  though  the 
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foliage.  We  make  use  of  the  measured  data  to  eompare  with  our  ray  traeing  model  in 
Chapter  IV. 
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III.  RAY  TRACING  MODEL 


In  this  chapter,  we  list  the  assumptions  and  considerations  of  the  ray  tracing 
model  and  discuss  the  formulation  of  the  model.  Next,  we  elaborate  on  the  modification 
the  model  suggested  in  [14]  for  computing  the  total  electric  field  for  the  GO  wave 
component.  A  set  of  general  equations  used  to  compute  the  path  length  for  each  GO  wave 
with  up  to  N  reflections  within  the  forest  layer  is  described.  A  lateral  wave  component 
introduced  by  [12]  is  described  and  examined  for  its  applicability  to  both  polarizations. 

A  three-layer  deterministic  model  is  generated  using  MATLAB  to  compute  the 
radiowave  propagation  loss.  This  model  is  setup  for  Configuration  2  from  Figure  3, 
where  both  the  transmitting  antenna  at  T  and  observation  point  located  at  O  are  immersed 
within  the  forest  layer  as  shown  in  Figure  6. 


Air  /^o’^o 


A  B 


Ground 

. ==> 

Reflected  waves 

- > 

Direct  wave 

- > 

Lateral  wave 


Figure  6.  Ray  tracing  model. 
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A,  ASSUMPTIONS  AND  CONSIDERATIONS 

For  the  computation  of  the  multi-bounce  waves  within  the  forest  layer,  reflections 
from  the  edges  of  forest  layer  are  negligible  and,  thus,  not  taken  into  consideration  [18]. 
This  assumption  is  valid  when  the  transmitting  antenna  and  observation  point  are  far 
away  from  the  edge  of  the  forested  environment. 

Only  the  primary  lateral  wave  contribution  and  the  single  and  double  bounce  rays 
are  taken  into  consideration  for  this  model.  As  discussed  previously,  the  longer  reflected 
path  distance  that  occurs  with  each  bounce  within  the  forest  attenuates  the  longer 
reflected  paths  more  than  the  shorter  reflected  paths.  Moreover,  with  the  multiplication  of 
additional  reflection  coefficient  factors,  we  can  assume  that  higher  orders  of  the  reflected 
waves  (more  than  one  bounce  off  the  same  boundary)  can  be  ignored. 

In  order  to  model  the  forest  as  a  uniform,  continuous  medium,  the  wavelength 
must  be  sufficiently  large  compared  to  the  forest  depth.  Wavelength  is  dependent  on  the 
frequency  of  operation,  and  the  maximum  usable  frequency  of  interest  for  our  model  is 
100  MHz,  which  is  the  same  criterion  as  Tamir  [12]. 

The  complex  permittivity  of  the  forest  layer  determines  the  refraction  angle  at  the 
boundary  layers  between  the  forest  and  air.  The  real  part  of  the  complex  permittivity  is 
used  to  determine  this  angle. 

There  is  potentially  an  ionospheric  wave  contribution,  but  it  is  not  taken  into 
consideration  in  the  computation  as  it  is  negligible  compared  to  the  lateral  waves.  This 
assumption  is  valid  when  the  distance  between  the  transmitter  and  observation  point  is 
less  than  100  km  [20]. 

B,  MODEL  FORMULATIONS 

The  ray  tracing  model  gives  the  total  field  strength  of  wave  components  emitted 
from  the  transmitting  antenna  at  the  observation  point.  The  computed  values  of  field 
strength  produced  by  the  propagation  mechanisms  described  in  Chapter  II  are  summed  at 
the  observation  point.  The  first  component  computes  the  field  strengths  from  both  direct 
and  the  multipath  rays  bounced  within  the  forest  layer,  and  the  second  component 


24 


computes  the  total  field  strength  generated  by  the  lateral  waves  formed  via  the  ground 
bounee  indueed  and  primary  lateral  wave  from  the  souree. 

A  set  of  equations  was  generated  to  compute  the  path  length  between  the 
transmitter  and  observation  point.  We  explain  the  modifieation  of  the  EM  wave  equation 
from  [14]  for  the  GO  component  in  this  section. 


1,  Propagation  Constant 


From  the  electromagnetic  theory  [27],  the  complex  propagation  constant  of  EM 
waves  for  a  lossy  medium  consists  of  real  and  imaginary  parts,  where  the  real  part  is  the 
attenuation  eonstant  and  the  imaginary  part  is  the  phase  eonstant.  Substituting  the 
eomplex  permittivity  of  the  forest  medium  Sj-  into  Eq.  (10),  we  find  that  the  real  part  of 


the  eomplex  propagation  eonstant,  that  is,  the  attenuation  eonstant  in  (Np/m)  of  the  forest 
layer,  ean  be  expressed  as  [26]; 


af  = 


1 


a 


1/Vao^o 


^fr  + 


f  (Jf  ^ 


) 


1/2 


a/2 


>-£ 


> 


(29) 


where  juq  and  Sq  are  the  permeability  and  permittivity  for  free  space,  respeetively.  The 


imaginary  part  is  the  phase  eonstant  j3 y  (rad/s)  of  the  forest  layer  is  expressed  as  [26]; 


CD 


1/Vao^O 


Sfr  + 


f  a 


nl/2 


V‘^^0  J 


sxn 


>  +  s 


fr 


(30) 


2.  Path  Length  Computation 

There  are  an  infinite  number  of  waves  that  bounce  off  the  forest-air  and  the  forest 
ground  boundary.  There  are  many  applieations  that  can  use  the  ray  traeing  method  to 
determine  effeets  in  the  communication  link  such  as  fading  and  multipath.  Our 
application  to  model  the  reflected  waves  within  the  forest  layer  is  less  complex  than 
many  other  multipath  problems  because  we  have  only  two  parallel  interfaces  to  consider. 
A  set  of  general  equations  is  derived  in  the  process  of  determining  the  path  length 
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between  the  transmitter  and  observation  point.  Let  m  be  the  number  of  bounces  off  the 
forest-air  boundary  and  n  be  the  number  of  bounces  off  the  ground-forest  boundary.  An 
example  of  four  possible  cases  in  computing  the  path  length  is  shown  in  Figure  7.  The 
allowed  ray  trajectories  are  obtained  by  forcing  Snell’s  law  to  be  satisfied  at  each 
reflection  point. 


Case  1  /»>/; 


Case  3  //;  =  /; 


Case  4  m  =  n 


Figure  7.  Four  possible  cases  for  the  path  length  computation. 


For  case  1  where  m>n  : 


^3  = 


D  =  xi  +  X2{m  +  n  -  \)  +  xt, 
H-hi  _  H  _H-h2 

xi  X2  X2 

D 

H-hl  .  “ 

- \-{m  +  n-l) - hi 

H-h2  ^  H-h2 


(31) 

(32) 

(33) 


(34) 
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For  case  2  where  m<n  : 


Z)  =  xj  +  X2  ( w  +  n  - 1)  +  X3 


H  ho 


^3  = 


Xj  X2  X3 

D 


— +  (m  +  n-l)  — +  1 

h2  ^  ’h2 


^fmn  -{m  +  n-\)^lH^  ■^X'l  +  ^  ^ ^2  +  -^3 

For  ease  3  where  m  =  n  : 

Z)  =  xj  +  X2  ( w  +  n  - 1)  +  X3 

ZZ-^i  _H  _h2 


Xj 


^3 


X2  X3 

D 


H-hl  ,  ..H  , 

- hlm  +  n-l) - hi 

/^2  ^  ^h2 


^fmn  =(w  +  n-l)'^ZZ^  +X2^  +"^(^“^1)  +^1^2^  +^3 

Finally,  for  ease  4  where  m  =  n  : 

Z)  =  xj  +  X2  ( w  +  n  - 1)  +  X3 


^3 


ZZ  ZZ  -  /?2 

Xi  X2  X3 

D 


hi  .  ,x  ZZ  , 

+  (m  +  n-lj - : — hi 


H-h2 


H-ho 


(35) 

(36) 

(37) 

(38) 

(39) 

(40) 

(41) 

(42) 

(43) 

(44) 

(45) 


^fmn ~hl)  +-^3  •  (46) 

These  four  sets  of  equations  are  used  to  eompute  the  refleeted  path  lengths 
between  the  transmitting  antenna  and  observation  point. 
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3. 


Geometrical  Optics  Waves  Formulation 


Electric  fields  generated  from  GO  waves  that  bounce  off  of  the  ground  plane  and 
forest- air  boundary  back  into  the  forest  towards  the  observation  point  are  defined  by  [14]: 


^GO  -  J 


JfTjflAz  e  (r^)'”(rg)” 


(47) 


mn 

where  is  wave  impedance  of  the  forest  layer,  /Az  is  the  source  dipole  moment  where 
/  is  the  current  flowing  through  the  dipole,  and  Az  is  the  length  of  dipole.  We  define  T y 
and  Eg  as  the  reflection  coefficients  for  the  forest-air  and  forest-ground  boundaries, 
respectively. 


4,  Antenna  Pattern  Factor 

The  derivation  of  Eq.  (47)  does  not  take  into  account  of  the  antenna  pattern  factor 
sin  9  illustrated  in  Eigure  8.  Overestimation  of  the  electric  field  can  occur  when  pattern 
factor  is  ignored.  This  overestimation  of  the  electric  field  is  only  applicable  for  the  V-pol 
for  a  vertically  oriented  dipole.  The  electric  field  for  H-pol  is  omnidirectional  and,  thus, 
is  independent  of  the  sin  9  far-field  pattern. 

The  far-field  distance  is  dependent  on  the  frequency  of  operation  and  the 
dimension  of  the  antenna  aperture.  In  order  to  determine  the  far-field  distance,  three 
conditions  must  be  satisfied.  The  expressions  for  the  three  criteria  are  [33]: 


n  2Az^ 

(48) 

R> 

A, 

R  >  5Az 

(49) 

R>\.6A 

(50) 

where  R  is  the  distance  from  dipole.  All  the  units  used  for  the  quantities  are  in  meters. 
Upon  computation  of  each  condition,  we  select  the  largest  distance  amongst  the  three 
values  as  Rg .  AX  distances  greater  than  Rjj-  from  transmitting  antenna,  the  EM  wave 

fronts  are  those  of  a  spherical  wave,  which  can  be  approximated  locally  as  a  plane  wave. 
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E-  plane  radiation 
pattern  for  £;i 


H-  plane  radiation 
pattern  for  l^il 


Figure  8.  E  -  plane  and  H  -  plane  radiation  plots  for  an  ideal  dipole  (after 

[33]). 


The  electric  field  in  free  space  that  is  generated  from  an  ideal  dipole  antenna  is 


[33]; 


^  .  /Az  ~ 

E  =  jcoju - Sint'  & 


(51) 


4;r  R 

where  R  is  the  distance  from  the  dipole.  The  relation  of  the  impedance  of  free  space  t/q 
to  the  phase  constant  /?  is  expressed  as: 

co/u 


% 


P 


(52) 


The  forest  layer  is  a  lossy  medium;  therefore,  we  replace  the  following  terms  co/^ , 
R,  and  -jp  in  Eq.  (51)  with  rjfpf,  and  -(^afRf„^„  +jj3fRfj„„Y  respectively. 

The  expression  for  computing  the  electric  field  for  the  GO  wave  is: 

-CCf  Rf  ~  j Pf  Rf  mn 


^mn  jp  f  P f 


lEz 


An 


Rfm 


sin  6*  . 


(53) 


The  difference  between  Eq.  (47)  and  Eq.  (53)  is  evident.  Firstly,  there  is  an 
assumption  made  in  Eq.  (47)  that  the  separation  distance  between  the  transmitting 
antenna  and  observation  point  is  very  large,  making  the  angle  9  close  to  90° .  Thus, 

sin  6*  Si  1 .  To  illustrate  the  difference,  we  examine  the  case  with  the  geometry  depicted  in 
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li  field  strength  (dlW-  m)  [•  Held  strength  (dBV'm) 


Figure  9.  For  simplification  of  the  problem,  we  compute  the  electric  fields  of  an  ideal 
dipole  in  free  space  near  a  ground  plane.  A  numerical  computation  based  on  Eq.  (47)  and 
Eq.  (53)  for  vertical  polarization  is  plotted  in  Figure  10. 


Air 


Reflected  waves 


Ground-air  boundary 
Ground 


Direct  wave 

Figure  9.  Geometry  setup  for  antenna  factor  analysis. 


Figure  10.  Effect  of  the  antenna  element  factor  on  direct  and  reflected  E  field  at 

100  MHz. 
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From  Figure  10,  the  impaet  of  the  antenna  pattern  on  the  eleetrie  field  at  near 
distanees  is  greater  for  the  refieeted  fields  eompared  to  the  direet. 


5.  Lateral  Wave  Formulation 

The  eomponents  of  the  lateral  wave  appear  in  the  form  [12]; 

-Icos^cos^+sin^j  (54) 

Ep=Ei  {^e -ijcos^cos^y +  -ijsin^  (55) 

E^  =  El  sin  ^  cos  y/  (56) 

where  y/  is  the  inclination  angle  of  the  dipole  antenna  from  horizontal  in  x  -  z  plane  and 
(j)  is  the  angle  from  the  x-z  plane.  The  term  E^  used  in  Eqs.  (54)  -  (56)  is  expressed 
as  [12]; 

-jpiD+sJsf  -l) 

60/Az  e  '  ^  ' 


El  = 


-1  D- 

where  s  is  the  total  path  length  in  the  foliage  between  the  transmitter  and  observation 
point.  This  path  length  is  approximated  as  s  -2H -hi  -/?2 .  The  complex  permittivity  of 
the  forest  layer  Sj-  is  defined  in  Eq.  (15).  As  we  are  interested  in  computing  the  electric 
field  of  a  wave  propagating  in  the  x-z  plane,  only  the  terms  E^  and  E^  are  required. 
The  E^  component  from  Eq.  (54)  is  used  for  computing  electric  field  for  V-pol 
=  90°,  ^  =  0°) ,  and  Ea  from  Eq.  (56)  is  used  in  computing  the  electric  field  for 


(57) 


H-pol  (^  =  0°,  ^  =  90°).  The  result  is  exactly  the  same  as  Eq.  (57);  therefore,  the 

formula  used  to  compute  the  lateral  wave  in  this  model  is  determined  by  Eq.  (57)  for  both 
polarizations. 
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6. 


Transmission  Loss  Definition 


The  computation  of  the  transmission  loss  TLj^  in  dB  at  a  distance  D  in  meters  is 
defined  as: 


7Z^(J5)  =  201ogio  Ej)  -201ogio 


meter 


(58) 


where 


D  and 


is  the  received  field  strength  at  the  observation  point  at  separation  distance 


^1  meter 


is  the  reference  field  strength  at  1 .0  m  from  transmitting  antenna  in  free 


space  without  ground,  as  shown  in  Figure  12. 


Transmitting  Observation  Point 

Antenna 

.  1  meter 

Both  height  of  the  antenna 
and  observation  point  are 

In  free  space 


Figure  1 1 .  Reference  electric  field  for  computing  the  transmission  loss. 

In  this  chapter,  we  presented  the  three-layered  model  and  a  homogeneous  model 
using  the  ray  tracing  method.  The  GO  component  and  the  lateral  wave  component  are 
considered  in  the  model.  Modification  of  the  GO  equation  was  made  to  improve  the  near 
distance  result  by  adding  the  antenna  factor  for  V-pol.  General  equations  were  derived  to 
compute  the  path  length  in  foliage  Rfmn  ■  The  definition  of  transmission  loss  was  given. 

Equation  (57)  is  used  to  compute  the  lateral  wave  for  both  polarizations.  In  Chapter  IV, 
the  ray  tracing  model  is  validated  by  using  the  results  generated  by  the  model  and  by 
comparing  those  results  with  measured  data. 
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IV.  RESULTS  AND  ANALYSIS 


In  this  chapter,  we  deseribe  the  two  eases  used  to  evaluate  the  ray  traeing  model. 
Results  are  diseussed  for  eaeh  test  ease. 

For  Case  I,  measured  data  extraeted  from  [14]  is  used.  The  simulation  setup  is 
done  sueh  that  the  height  of  the  transmitting  antenna  remains  at  3.96  m  for  all 
measurements,  and  the  height  of  the  observation  point  varies  from  5.0  m  to 
28.96  m.  Measured  data  were  taken  for  both  V-pol  and  H-pol  at  frequeneies  of 
25  MHz,  50  MHz,  and  100  MHz. 

For  Case  II,  we  validate  the  ray  traeing  model  using  two  empirieal  models 
mentioned  in  Chapter  II  (the  JB  model  and  Tewari’s  empirieal  model)  together  with 
measured  data  extraeted  from  [7].  The  50  MHz  data  point  from  [7]  is  eompared  with  the 
models.  The  root  mean  square  (RMS)  error  in  dB  introdueed  in  [31]  is  used  to 

define  a  fitting  value  to  quantify  the  performanee  of  the  ray  traeing  model  with  respeet  to 
measurements. 

A.  CASE  I  SETUP  AND  SIMULATION 

Two  sets  of  parameters  are  defined  for  this  test  setup.  First  is  the  geometry  for  the 
setup,  whieh  ineludes  the  physieal  heights  of  the  forest,  transmitting  antenna,  and 
observation  point,  and  the  separation  distanee.  The  geometry  inputs  for  the  setup  are 
listed  in  Table  5. 


Table  5.  Geometrieal  inputs  for  Case  I. 


Parameter 

Value 

Units 

Average  Height  of  trees  in  forest  layer,  H 

30.48 

m 

Separation  Distanee,  D 

1600 

m 

Height  of  transmitting  antenna,  h\ 

3.96 

m 

For  H-pol 

Height  of  observation  point,  ^2 

5-29 

m 

For  V-pol 

Height  of  observation  point,  ^2 

7-29 

m 
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The  next  set  of  input  parameters  are  the  eleetrieal  properties  of  eaeh  medium.  The 
ground  relative  permittivity  and  ground  eonduetivity  listed  in  Table  6  is  with  referenee  to 
[14],  These  values  are  typieal  ground  values  as  reported  in  [20].  For  the  eleetrieal 
properties  of  the  forest  media,  we  refer  to  Eq.  (15)  from  Chapter  II,  where  the  forest 
layer’s  effeetive  permittivity  and  eonduetivity  cry  ean  be  eombined  to  form  a 

eomplex  permittivity. 

The  forest  layer  is  modeled  as  an  anisotropie  medium  in  [14].  It  has  an  effeetive 
permittivity  and  eonduetivity  for  the  z-direetion  as  well  as  the  tangential  direetion.  The 
eleetrieal  properties  of  the  z-direetion  primarily  affeets  the  V-pol  waves  propagating  in 
the  forest  layer.  The  eleetrieal  properties  for  the  tangential  direetion  primarily  affeet  the 
H-pol  waves  propagating  in  the  forest  layer.  The  ray  traeing  model  uses  homogeneous 
media  for  all  layers;  thus,  only  the  eleetrieal  properties  eorresponding  to  the  polarization 
of  interest  are  used  as  input.  So  for  V-pol,  the  effeetive  relative  permittivity  is  represented 
by  and  eonduetivity  by  cr^ .  The  effeetive  relative  permittivity  and  eonduetivity  for 

H-pol  is  represented  by  and  cr^ ,  respeetively.  The  initial  values  used  in  the 
simulation  are  the  ones  suggested  in  [14]  and  are  listed  in  Table  6. 


Table  6.  Initial  eleetrieal  properties  of  media  for  Case  I. 


Parameter 

Value 

Units 

Frequeney 

25,  50,  100 

MHz 

Relative  permittivity  of  air  medium. 

1 

- 

eonduetivity  of  air  medium 

0 

mS/m 

Relative  permittivity  of  ground  medium 

15 

- 

eonduetivity  of  ground  medium 

10 

mS/m 

Relative  permittivity  of  forest  medium  (V-pol) 

1.0530 

- 

eonduetivity  of  forest  medium  (V-pol) 

0.1180 

mS/m 

Relative  permittivity  of  forest  medium  (H-pol) 

1.0080 

- 

eonduetivity  of  forest  medium  (H-pol) 

0.0300 

mS/m 

Upon  establishing  all  the  neeessary  values  for  Case  I,  an  iterative  proeess  is 

performed  to  obtain  the  modified  effeetive  parameters  that  give  the  best  agreement 

between  the  simulated  and  measured  data.  The  proeess  eonsists  of  two  steps.  Firstly,  the 
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relative  permittivity  of  the  forest  layer  is  fixed  at  a  eonstant  value  while  ehanging  the 
value  of  eonduetivity  of  the  forest  media.  The  initial  eonstant  value  of  relative 
permittivity  for  the  forest  layer  used  is  with  referenee  to  optimized  data  in  [14],  Upon 
getting  a  minimum  average  error  value,  the  relative  permittivity  value  is  then  tuned  to 
aehieve  lowest  possible  average  error. 

1.  Case  I  Results 

The  final  adjusted  relative  effeetive  permittivity  and  eonduetivity  cr^ 

eorresponding  to  the  z-direetion  are  tabulated  in  Table  7  at  three  frequeney  points 
(25  MHz,  50  MHz,  and  100  MHz).  Similarly  in  Table  8,  the  final  adjusted  relative 
effeetive  permittivity  and  eonduetivity  cjf  that  eorrespond  to  the  tangential  direetion 

are  tabulated. 

The  simulated  loss  versus  height  eurves  for  frequeneies  of  25  MHz,  50  MHz,  and 
100  MHz  for  V-pol  are  plotted  in  Figure  12  to  Figure  14,  respeetively.  Similarly,  the 
simulated  loss  eurves  for  the  same  frequeney  points  for  H-pol  are  plotted  in  Figure  15  to 
Figure  17,  respeetively.  In  both  oases,  measured  data  points  are  superimposed  on  the 
simulated  results. 


Table  7.  Effeetive  permittivity  and  eonduetivity  for  the  ray  traoing  method 

and  Sommerfeld  solution  from  [14]  (V-pol). 


V-pol 

Sommerfeld  solution 

Ray  tracing  model 

Frequeney  (MHz) 

(mS/m) 

(mS/m) 

25 

1.0530 

0.1180 

1.0600 

0.1010 

50 

1.0180 

0.0730 

1.0400 

0.0930 

100 

1.0060 

0.0420 

1.0250 

0.0760 
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Table  8.  Effective  permittivity  and  conductivity  for  the  ray  tracing  method 

and  Sommerfeld  solution  from  [14]  (H-pol). 


H-pol 

Sommerfeld  solution 

Ray  tracing  model 

Frequency  (MHz) 

^4 

(mS/m) 

^4 

(mS/m) 

25 

1.0080 

0.0300 

1.0090 

0.0300 

50 

1.0100 

0.0370 

1.0110 

0.0410 

100 

1.0100 

0.0420 

1.0610 

0.0550 

Height  of  Observation  Point,  (m) 


Figure  12.  Measured  data  from  [14]  versus  ray  tracing  model  (1.6  km,  25  MHz, 

V-pol). 
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Height  of  Observation  Point,  (m) 


Figure  13.  Measured  data  from  [14]  versus  ray  tracing  model  (1.6  km,  50  MHz, 

V-pol). 


Figure  14.  Measured  data  from  [14]  versus  ray  tracing  model  (1.6  km,  100 

MHz,  V-pol). 
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Figure  15.  Measured  data  from  [14]  versus  ray  tracing  model  (1.6  km,  25  MHz, 

H-pol). 


Figure  16.  Measured  data  from  [14]  versus  ray  tracing  model  (1.6  km,  50  MHz, 

H-pol). 
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Figure  17.  Measured  data  from  [14]  versus  ray  tracing  model  (1.6  km,  100 

MHz,  H-pol). 


2,  Case  I  Analysis  and  Comments 

The  results  generated  by  the  ray  tracing  model  were  based  on  the  assumption  that 
forest  media  is  isotropic  and  homogenous.  The  average  difference  between  each 
measured  data  point  with  the  simulated  data  point  at  each  frequency  is  tabulated  in 
Table  9.  The  difference  for  both  polarizations  was  within  ±  1  dB,  similar  to  that  of  the 
Sommerfeld  solution  [14].  The  idea  of  using  a  homogenous  medium  to  represent  the 
forest  layer  provided  a  satisfactory  correlation  to  the  measured  data  (within  1  dB)  up  to  a 
frequency  of  100  MHz.  Another  observation  was  the  reduction  of  loss  with  increase  in 
height  of  observation  point.  This  is  expected  because  the  path  length  in  the  lossy  foliage 
decreases  as  the  antenna  height  increases. 
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Table  9.  Average  error  between  measured  data  from  [14]  and  ray  tracing 

model  with  adjusted  parameter. 


Average  error  compared  between  model  and  measured  data  (dB) 

H-pol  to  H-pol 

V-pol  to  V-pol 

Frequency 

(MHz) 

Sommerfeld 

solution 

Ray  Tracing 
Model 

Sommerfeld 

solution 

Ray  Tracing 
Model 

25 

0.60 

0.82 

0.30 

0.90 

50 

0.80 

0.23 

0.40 

0.48 

100 

0.90 

0.97 

0.40 

0.72 

Up  to  this  point,  the  model  was  validated  to  100  MHz  at  a  separation  distance  of 


1.6  km.  Further  testing  is  needed  to  validate  the  range  applicability  of  the  ray  tracing 
model. 


B.  CASE  II  SETUP  AND  SIMULATION 

In  Case  I,  we  commented  that  the  ray  tracing  model  was  able  to  provide 
satisfactory  estimation  of  the  transmission  loss  at  a  fixed  separation  distance  of  up  to 
100  MHz.  Further  verification  on  the  applicability  of  the  ray  tracing  model  was  made  by 
comparing  it  with  another  set  of  measured  data  collected  for  the  northern  part  of  India 
from  [7].  Procedures  established  in  Case  I  were  repeated,  where  we  defined  the  geometry 
inputs  (refer  to  Table  10)  and  the  estimated  electrical  properties  of  each  medium  (refer  to 
Table  11).  For  Case  II,  the  electrical  properties  of  the  forest  were  taken  from  the 
measured  data  via  the  inverse  method  as  documented  in  [26]. 


Table  10.  Geometrical  inputs  for  Case  II. 


Parameter 

Value 

Units 

Average  Height  of  trees  in  forest  layer,  H 

20 

m 

Separation  Distance,  D 

1-4200 

m 

Height  of  transmitting  antenna, 

5 

m 

Height  of  observation  point,  /i2 
(For  H-pol  or  V-pol) 

5 

m 
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Table  1 1 .  Electrical  properties  of  each  medium  for  Case  II. 


Parameter 

Value 

Units 

Frequency 

50 

MHz 

Relative  permittivity  of  air  medium, 

1 

- 

Conductivity  of  air  medium 

0 

mS/m 

Relative  permittivity  of  ground  medium 

15 

- 

Conductivity  of  ground  medium 

10 

mS/m 

Relative  permittivity  of  forest  medium 

1.065 

- 

Conductivity  of  forest  medium 

0.135 

mS/m 

1.  Case  II  Results 

In  conjunction  with  the  measured  data  used  [7],  two  empirical  models  were  used 
for  comparison.  The  data  are  plotted  in  Figure  18  and  Figure  19. 


V-pol. 
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Figure  19.  Measured  data  from  [7]  versus  prediction  from  models  at  50  MHz, 

H-pol. 


A  useful  quantity  to  compare  the  performance  of  the  models  is  needed.  A  good 
figure  of  merit  is  the  RMS  error  defined  in  [31]; 


^RMS 


(59) 


where  is  the  RMS  error  in  dB,  N  is  number  of  measured  data  points,  and  £]•  is  the 


absolute  difference  in  dB  between  the  measured  and  the  simulated  data.  The  RMS  errors 
for  both  polarizations  at  50  MHz  are  tabulated  in  Table  12. 


Table  12.  RMS  error  at  50  MHz. 


Polarization 

Ray  Tracing  Model, 
^RMS  (dB) 

Tewari’s  Empirical 
Model,  (dB) 

JB  Empirical 
Model,  Srms  (dB) 

V-pol 

3.7685 

1.1873 

1.2047 

H-pol 

2.6869 

0.6528 

0.7097 
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2,  Case  II  Analysis  and  Comments 

The  RMS  error  value  of  the  ray  tracing  model  was  much  higher  compared  to  the 
other  two  empirical  models.  From  [12],  it  was  suggested  that  the  formulation  for  the 
transmission  loss  only  applies  for  separation  distance  D  greater  than  1.0  km.  To 
investigate  this,  we  considered  only  measured  data  points  where  separation  distance  D 
was  larger  than  1 .0  km.  A  new  computation  of  the  RMS  error  was  performed  and  the 
results  tabulated  in  Table  13. 


Table  13.  RMS  error  at  50  MHz  with  D  greater  than  1.0  km. 


Polarization 

Ray  Tracing  Model, 
^RMS  (dB) 

Tewari’s  Empirical 
Model, 

^RMS  (dB) 

JB  Empirical 
Model, 

^RMS  (dB) 

V-pol 

2.7073 

2.4377 

2.5979 

H-pol 

4.3222 

1.3141 

1.3813 

RMS  error  for  V-pol  is  reduced  and  closer  to  both  empirical  models  when  the 
simulated  data  are  compared  to  measured  data  at  separation  distance  D  greater  than 
1 .0  km;  however,  we  still  observe  large  differences  for  H-pol.  This  could  be  due  to  two 
possible  reasons.  First,  in  the  formulation  of  the  lateral  wave  equation,  both  the 
polarizations  use  the  same  expressions.  Second,  the  electrical  properties  for  both 
V-pol  and  H-pol  are  of  the  same  values  with  reference  to  Table  2.  Based  on  findings  in 
[14],  a  higher  forest  conductivity  cry  value  for  V-pol  is  used  compared  to  the  H-pol.  This 

is  explained  by  stronger  scattering  of  the  vertical  polarized  electric  fields  due  to  the 
vertical  tree  trunks  in  the  forest  layer  [14].  To  further  investigate  this  characteristic,  the 
forest  conductivity  was  reduced  by  30  percent  from  0.135  mS/m  to  0.095  mS/m.  The  new 
resultant  RMS  error  for  the  ray  tracing  model  is  tabulated  in  Table  14. 
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Table  14.  RMS  error  with  D  greater  than  1.0  km  and  lower  forest 

eonduetivity. 


Polarization 

Ray  Tracing  Model, 
^RMS  (dB) 

Tewari’s  Empirical 
Model, 

^RMS  (dB) 

JB  Empirical 
Model, 

^RMS  (dB) 

H-pol 

1.3415 

1.3141 

1.3813 

The  new  resultant  RMS  error  based  on  the  new  forest  eonduetivity  is  redueed. 
This  eonstitutes  a  smaller  differenee  between  the  other  two  empirieal  models  at  a 
separation  distanee  greater  than  1.0  km.  This  finding  suggests  that  higher  attenuation  is 
experieneed  by  the  vertieal  polarized  eleetrie  field  eompared  to  horizontal  polarized 
eleetrie  fields  when  propagating  within  the  forest  media;  however,  this  reduction  in  forest 
conductivity  should  not  be  taken  as  a  proportional  effect  compared  to  the  scattering 
effects  from  the  tree  trunks.  It  only  means  that  for  estimation  of  transmission  loss  within 
a  forest  media,  the  attenuation  by  the  V-pol  electric  field  is  in  general  greater  than  that  of 
the  H-pol  electric  field. 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 


A,  SUMMARY 

A  ray  tracing  model,  based  on  the  three-layer  planar  stratified  medium  shown  in 
Figure  2,  was  coded  in  MATLAB.  It  allows  users  to  numerically  compute  the  electric 
field  and  to  analyze  direct  waves,  multiple  refleeted  waves,  and  lateral  waves  or  a 
summation  of  all  of  these  eomponents 

For  the  formulation  to  compute  the  electric  field  for  refleeted  and  direct  rays,  an 
antenna  pattern  faetor  was  introduced  [14],  It  was  found  that  V-pol  reflected  waves  are 
affeeted  more  by  the  faetor  sin  9  eompared  to  the  direet  waves.  The  antenna  pattern 
faetor  is  not  applieable  to  the  H-pol  eleetrie  fields  as  the  antenna  pattern  is 
omnidireetional.  A  set  of  equations  was  introduced  to  compute  the  multi-bounee  path 
lengths  within  the  foliage.  These  equations  were  derived  by  applying  the  well-known 
Snell’s  law  of  reflection  [27]  as  well  as  the  properties  of  symmetrieal  triangles  where  the 
all  sides  of  eaeh  triangle  are  proportional.  For  the  eomputation  of  the  lateral  waves,  the 
equation  derived  in  [12]  was  used  and  is  applicable  for  both  V-pol  and  H-pol. 

The  ray  traeing  model  was  eompared  to  measured  data  for  two  cases.  In  Case  I, 
simulated  results  at  25  MHz,  50  MHz,  and  100  MHz  for  both  H-pol  and  V-pol  were 
compared  to  the  measured  results  reported  in  [14].  Results  suggested  that  the 
homogeneous  media  representing  the  forest  layer  was  applicable  up  to  100  MHz  and  at  a 
separation  distance  of  1.6  km.  This  upper  frequency  limit  of  the  ray  tracing  model  was 
consistent  with  the  applicable  frequency  reported  in  [12].  In  Case  II,  another  comparison 
was  performed  to  determine  the  useful  range  of  the  model.  A  comparison  at  50  MHz  was 
made  with  two  other  empirical  models  (JB  model  and  Tewari’s  empirical  model)  together 
with  measured  data  approximated  from  [7].  To  quantify  the  performance,  a  RMS  error 
defined  in  [31]  was  introduced.  The  RMS  error  of  the  ray  tracing  model  for  V-pol  at 
50  MHz  showed  good  agreement  with  two  other  emp ideal  models  (the  JB  model  and  the 
Tewari’s  empirieal  model)  at  separation  distanees  beyond  1.0  km.  For  the  H-pol  at  the 
same  frequeney,  we  observed  large  differences  in  the  RMS  error  between  the  ray  traeing 
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model  and  the  other  two  empirical  models.  This  was  explained  by  the  scattering  effect  of 
the  tree  trunks.  It  was  found  that  the  V-pol  electric  fields  experienced  more  attenuation 
compared  to  the  H-pol  electric  fields  due  to  the  greater  scattering  effect  caused  by  the 
physical  vertical  tree  trunks  within  the  forest  [14].  This  was  confirmed  by  tuning  the 
conductivity  of  the  forest  media  to  a  lower  value  that  corresponds  to  lower  transmission 
loss. 

B,  CONCLUSION 

In  conclusion,  a  physics-based  model  provides  a  fundamental  knowledge  of  wave 
propagation  in  stratified  media  useful  for  further  research  purposes.  A  propagation  loss 
prediction  model  using  the  ray  tracing  method  was  successfully  implemented  for 
Configuration  2  in  Figure  3.  Leveraging  on  the  previous  works  done  by  Tamir  [18]  and 
Li  [14],  we  were  able  to  formulate  the  model  for  applications  where  communication 
devices  are  operating  at  a  frequency  less  than  100  MHz,  with  the  minimum  distance 
between  the  communication  devices  being  at  least  1.0  km  apart  within  a  forest 
environment.  This  ray  tracing  model  used  the  three-layer  planar  stratified  medium  shown 
in  Figure  2  and  was  coded  using  MATLAB  software.  It  allows  users  to  compute  the 
electric  field  for  direct  waves,  multiple  reflected  waves,  and  lateral  waves  or  a  summation 
of  all  these  waves  within  the  forest  medium  numerically.  Through  the  use  of  analytical 
methods,  the  geometry  of  the  problem  can  be  extended  to  greater  depth  of  the  forest 
compared  to  the  CEM  used  in  previous  work  [13],  where  the  computational  resources  and 
time  needed  are  greater. 

In  Case  II,  it  was  suggested  that  vertical  electric  fields  propagating  within  the 
forest  medium  in  general  experience  more  scattering  effects  due  to  the  vertical  tree 
trunks;  therefore,  higher  attenuations  occur  compared  to  the  horizontal  polarized  electric 
fields  [14]. 
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C.  RECOMMENDATIONS 

During  the  course  of  analyzing  the  simulated  transmission  loss  in  Case  I,  the 
effective  parameters  of  the  forest  layer  used  from  [14]  played  an  important  role.  These 
values  served  as  an  initial  input  parameter  and  were  subsequently  tuned  during  the 
iterative  process  to  match  the  measured  data.  There  is  a  need  to  determine  these  effective 
parameters  either  via  empirical  models  or  via  measurements  when  no  initial  input 
parameters  are  provided. 

The  work  done  in  this  thesis  serves  as  an  initial  platform  that  allows  further 
modification  of  the  formulation  to  support  other  types  of  configurations  depicted  in 
Figure  3.  This  is  essentially  useful  for  different  types  of  communication  systems 
deployed  in  various  scenarios  that  can  be  similar  to  the  possible  configurations 
mentioned.  The  formulation  of  the  model  could  reference  Tamir’s  ray  tracing  method  in 
[18]. 

Another  limitation  of  the  model  that  requires  further  improvement  is  short  range 
transmission  loss  prediction.  From  Case  I  and  Case  II,  the  model  is  applicable  in 
predicting  transmission  loss  for  communication  links  greater  than  1  km.  One  possible 
software  tool  for  modeling  of  the  short  range  propagation  loss  at  VHF  band  is  CST 
Microwave  Studio.  This  overcomes  the  limitation  of  the  existing  model;  but  the 
simulation  time  can  be  extremely  long  (days). 
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APPENDIX  A.  SOFTWARE  DESCRIPTION 


The  ray  tracing  model  is  coded  using  the  MATLAB  program  and  inputs  are  keyed 
into  the  script  window  before  running  the  codes.  There  are  several  inputs  required  from 
the  user,  such  as  the  dimensions  of  the  forest  in  terms  of  depth  and  average  height  of  the 
forest  layer.  The  list  of  parameters  is  populated  in  Table  15. 


Table  15.  Input  parameters  for  the  ray  tracing  model. 


Description 

Symbol  used  in 
MATLAB  code 

Units 

Average  Height  of  the  forest  layer 

H 

m 

Height  of  transmitting  antenna 

hi 

m 

Height  of  observation  point 

hi 

m 

Separation  Distance 

D 

m 

Frequency 

flvlHz 

MHz 

Polarization  (V-pol,  H-pol) 

POEARSel 

- 

Effective  relative  permittivity  of  forest  layer 

ef 

- 

Effective  relative  permittivity  of  ground 
layer 

eg 

- 

Effective  relative  permittivity  of  air  layer 

ea 

- 

Conductivity  of  forest  layer 

sigma  1 

S/m 

Conductivity  of  ground  layer 

sigma2 

S/m 

Conductivity  of  air  layer 

sigmaS 

S/m 

Current  flowing  through  dipole 

I 

A 

Eength  of  dipole 

deltaZ 

m 
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The  overview  of  the  ray  traeing  model  eoded  in  MATLAB  is  illustrated  in  Figure  20. 


MATLAB  simulation  environment 


Ray  tracing  model  coded  in  MATLAB  Script 

Manual  Input 
parameters 
from  Table  1  [■ 

into 

MATLAB 
Script 


> 


BLOCK  A 


BLOCK  B 


OUTPUT 

•  Transmission  Loss 

•  GO  wave  E  field 

•  Lateral  wave  E  field 


Figure  20.  Overview  of  ray  traeing  model  in  MATLAB  simulation  environment 


Figure  21 .  Generie  flow  diagram  within  Bloek  A. 
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BLOCK  B 


Figure  22.  Generic  flow  diagram  within  Block  B. 


The  ray  tracing  model  was  coded  in  MATLAB  that  allows  user  to  decide  and 
code  their  preferred  output  for  analysis.  In  this  thesis,  computation  of  transmission  loss 
over  different  observation  height  was  coded  for  Case  I.  For  Case  II,  the  transmission  loss 
over  separation  distance  was  computed. 
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APPENDIX  B.  ABSOLUTE  DIFFERENCE  DATA  FOR  CASE  I 


Table  16.  Absolute  difference  (dB)  between  ray  tracing  model  and  measured 

data  at  25  MHz  (V-pol). 


Height  of 
Observation  Point, 

h  (m) 

TLj^  (dB)  at 

25  MHz  from  [14] 

Ray  Tracing 
Model  (dB) 

Absolute  difference 
(dB) 

5 

-127.5 

-125.8 

1.7 

10 

-124.2 

-124.0 

0.2 

15 

-119.4 

-121.0 

1.6 

20 

-117.8 

-118.0 

0.2 

25 

-113.6 

-115.1 

1.5 

28.96 

-112.5 

-112.7 

0.2 

Table  17.  Absolute  difference  (dB)  between  ray  tracing  model  and  measured 

data  at  50  MHz  (V-pol). 


Height  of 
Observation  Point, 

^  (m) 

TLj:,  (dB)  at 

50  MHz  from  [14] 

Ray  Tracing 
Model  (dB) 

Absolute  difference 
(dB) 

5 

-132.2 

-132.4 

0.2 

10 

-128.6 

-130.2 

1.6 

15 

-126.7 

.126.7 

0.0 

20 

-123.1 

-123.1 

0.0 

25 

-120.0 

-119.6 

0.4 

28.96 

-117.5 

-116.8 

0.7 

Table  18.  Absolute  difference  (dB)  between  ray  tracing  model  and  measured 

data  at  100  MHz  (V-pol). 


Height  of 
Observation  Point, 

h  (m) 

TLj)  (dB)  at 

100  MHz  from  [14] 

Ray  Tracing 
Model  (dB) 

Absolute  difference 
(dB) 

5 

-137.5 

-135.7 

1.8 

10 

-133.3 

-133.4 

0.1 

15 

-129.4 

-129.6 

0.2 

20 

-126.9 

-125.8 

1.1 

25 

-121.9 

-122.0 

0.1 

28.96 

-120.0 

-119.0 

1.0 
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Table  19.  Absolute  difference  (dB)  between  ray  tracing  model  and  measured 

data  at  25  MHz  (H-pol). 


Height  of 
Observation  Point, 

h  (m) 

TLj^  (dB)  at 

25  MHz  from  [14] 

Ray  Tracing 
Model  (dB) 

Absolute  difference 
(dB) 

5 

-106.7 

-103.7 

3.0 

10 

-101.9 

-101.9 

0.0 

15 

-99.4 

-100.0 

0.6 

20 

-97.2 

-98.1 

0.9 

25 

-96.1 

-96.2 

0.1 

28.96 

-95.0 

-94.7 

0.3 

Table  20.  Absolute  difference  (dB)  between  ray  tracing  model  and  measured 

data  at  50  MHz  (H-pol). 


Height  of 
Observation  Point, 

h  (m) 

TLj:,  (dB)  at 

50  MHz  from  [14] 

Ray  Tracing 
Model  (dB) 

Absolute  difference 
(dB) 

5 

-116.7 

-116.7 

0.0 

10 

-112.5 

-113.9 

1.4 

15 

-111.1 

-111.1 

0.0 

20 

-108.3 

-108.3 

0.0 

25 

-105.6 

-105.6 

0.0 

28.96 

-103.3 

-103.3 

0.0 

Table  2 1 .  Absolute  difference  (dB)  between  ray  tracing  model  and  measured 

data  at  100  MHz  (H-pol). 


Height  of 
Observation  Point, 

h  (m) 

TLj^  (dB)  at 

100  MHz  from  [14] 

Ray  Tracing 
Model  (dB) 

Absolute  difference 
(dB) 

5 

-122.8 

-12.5 

0.7 

10 

-121.6 

-121.6 

0.0 

15 

-121.9 

-119.8 

2.1 

20 

-119.4 

-118.0 

1.4 

25 

-116.1 

-116.1 

0.0 

28.96 

-113.1 

-114.7 

1.6 
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APPENDIX  C.  ABSOLUTE  DIFFERENCE  DATA  FOR  CASE  II 


Table  22.  Absolute  difference  (dB)  between  ray  tracing  model  and  measured 

data  at  50  MHz  (V-pol). 


Separation 
Distance,  D  (m) 

TLj^  (dB)  at 

50  MHz  from  [7] 

Ray  Tracing 
Model  (dB) 

Absolute 

difference, 

Ei  (dB) 

40 

-63.2 

-40.4 

22.8 

90 

-80.0 

-56.1 

23.9 

192 

-94.7 

-84.8 

9.9 

400 

-102.6 

-101.2 

1.4 

600 

-114.7 

-108.2 

6.5 

800 

-121.6 

-113.2 

8.4 

1000 

-127.4 

-117.1 

10.3 

2100 

-129.5 

-130.0 

0.5 

3100 

-134.7 

-136.8 

2.1 

4100 

-144.2 

-141.6 

2.6 

Table  23.  Absolute  difference  (dB)  between  ray  tracing  model  and  measured 

data  at  50  MHz  (H-pol). 


Separation 
Distance,  D  (m) 

TLj^  (dB)  at 

50  MHz  from  [7] 

Ray  Tracing 
Model  (dB) 

Absolute 

difference, 

Ei  (dB) 

40 

-54.7 

-40.3 

14.4 

91 

-65.3 

-63.9 

1.4 

192 

-80.0 

-82.2 

2.2 

400 

-92.6 

-101.0 

8.4 

600 

-102.1 

-108.2 

6.1 

800 

-103.2 

-113.2 

10.0 

1000 

-110.5 

-117.1 

6.6 

2100 

-118.9 

-130.0 

11.1 

3100 

-126.3 

-136.8 

10.5 

4100 

-136.8 

-141.6 

4.8 
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